Tracing masses of ground-state light-quark mesons 
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We describe a symmetry-preserving calculation of the light-quark meson spectrum, which com- 
bines a description of pion properties with reasonable estimates of the masses of heavier mesons, 
including axial-vector states. The kernels used in formulating the problem are essentially non- 
perturbative, incorporating effects of dynamical chiral symmetry breaking (DCSB) that were not 
previously possible to express. Our analysis clarifies a causal connection between DCSB and the 
splitting between vector and axial-vector mesons, and exposes a key role played by the anomalous 
chromomagnetic moment of dressed-quarks in forming the spectrum. 
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I. PROLOGUE 

Spectroscopy is basic to exposing the character of Na- 
ture's fundamental forces, e.g.: the hydrogen atom spec- 
trum is intimately connected with the development of 
quantum mechanics and quantum electrodynamics; and 
cataloguing the surfeit of hadrons - the pion- and proton- 
like particles, initially supposed to be elementary - led to 
the notion of quarks, with their fractional electric charge, 
and ultimately to quantum chromodynamics (QCD), the 
strongly-interacting piece of the Standard Model. 

QCD is peculiar owing to the empirical fact of con- 
finement] namely, whilst it is formulated in terms of 
quarks and gluons, which carry the color quantum num- 
ber, these degrees of freedom have never been directly 
observed. Determination of the spectrum of color-singlct 
hadrons is therefore one of the few ways by which to elu- 
cidate the long-range behavior of the interaction between 
quarks and gluons, and thereby grasp the essence of Na- 
ture's sole confining force. It might thus be surprising 
that more than thirty years after the discovery of the 
first light-quark meson with mass greater than 1 GeV, 
there is no reliable, symmetry-preserving computation of 
the ground-state spectrum of light-quark mesons. This 
is particularly disturbing given the importance of spec- 
troscopy at next-generation hadron physics facilities 

Confinement is much misapprehended. We therefore 
remark that the static potential measured in lattice-QCD 
simulations is not related in any known way to the puz- 
zle of light-quark confinement. Light-quark creation and 
annihilation effects are fundamentally nonperturbative. 
It is thus impossible in principle to compute a poten- 
tial between light quarks. Alternatively, arguments re- 
lating confinement to the analytic properties of QCD's 
Schwinger functions have been presented 0, from 
which perspective the question of light-quark confine- 
ment may be translated into the challenge of charting the 
infrared behavior of QCD's universal /3-function. Solv- 
ing this problem is a basic goal of modern physics, which 
can be addressed in any framework that enables the non- 
perturbative evaluation of rcnormalization constants. 



The light-hadron spectrum has challenged lattice-QCD 
since its formulation almost forty years ago The ap- 
proach involves a range of systematic errors, amongst 
which explicit chiral symmetry breaking is particularly 
damaging to studies of light-quark systems. It has only 
recently become possible to tame these errors and repro- 
duce the mass and width of the p- meson However, in 
connection with more massive light-quark mesons,avail- 
able results are much affected by lattice artifacts @, 0] ■ 

We prefer a continuum nonperturbative approach and, 
exploiting a novel form of the Bethc-Salpeter equation 
(BSE) [8[, formulate a tractable, symmetry preserving 
expression of the meson bound-state problem. This treat- 
ment connects the /3-function to observables, so that com- 
parison between computations and the spectrum may be 
used to constrain the /3-function's long-range behavior, 
as illustrated in Ref. H- 



II. BOUND-STATE EQUATIONS 

Our approach begins with the gap equation 

,A 

Sf{p)- 1 = Z 2 {i 1 -p + m) m ) + Z x \ g 2 D^{p - q) 

J da 
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where: / labels quark flavor; D^ v is the gluon propa- 
gator; T[, the quark-gluon vertex; J d , a Poincare in- 
variant rcgularization of the integral, with A the regu- 
larization scale; m bm (A), the current-quark bare mass; 
and Z 1j2 (CjA), respectively, the vertex and quark wave 
function renormalization constants, with £ the renormal- 
ization point - dependence upon which we do not usually 
make explicit. The solution is the quark propagator, 

Sfip)- 1 =i-f P A f (p 2 ) + B f (p 2 ). (2) 

It is obtained from Eq. ([1]) augmented by a renormaliza- 
tion condition. 

The mass of all mesons with the same quantum num- 
bers may be obtained from a single inhomogencous BSE. 
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The dressed-quark propagator is crucial in construct- 
ing its kernel. The gap equation involves the dressed- 
quark gluon vertex, T^. Following Ref. 0, one can now 
construct a symmetry-preserving kernel for the Bethe- 
Salpeter kernel given any form for T M . Hence, in solving 
bound-state problems one is no longer reliant upon step- 
by-step improvements in the computation of T M [10h12| ; 
instead, one may use all available, reliable information 
to construct the best possible Ansatz. This enables one 
to incorporate crucial nonpcrturbative effects, which any 
finite sum of contributions is incapable of capturing. In 
this way the nonpcrturbative phenomenon of dynamical 
chiral symmetry breaking (DCSB) was shown to generate 
material, momentum-dependent anomalous chromo- and 
electro-magnetic moments for dressed light-quarks (l3| . 

DCSB is a keystone of the Standard Model. This re- 
markable mass-generating mechanism can be explained 
via S(p). In the chiral limit, the mass function, M(p 2 ) = 
B(p )/A(p ), is identically zero at any finite order in per- 
turbation theory. However, DCSB generates mass from 
nothing. Thus, in chiral-QCD, despite the absence of 
an explicit mass source, M(p 2 ) is nonzero and stro ngly 
momentum-dependent, with M(0) w 0.5 GeV [HI. fl5f. 
DCSB is responsible for constituent-quark masses and in- 
timately connected with confinement but it is not known 
whether the connection is accidental or causal. Never- 
theless, DCSB is the most important mass-generating 
mechanism for visible matter in the Universe. It gen- 
erates roughly 98% of a proton's mass and amplifies the 
Higgs mechanism for explicit symmetry breaking (l6| . 

Herein we calculate masses of ground-state spin-zero 
and -one light-quark mesons, in order to illuminate the 
impact of DCSB on the spectrum. Given its central role, 
we explicate the problem in the pseudoscalar, J PC — 
h , and axial- vector, 1 ++ , channels and simply indicate 
the changes required to study ++ , l -- , l + ~. 

Pseudoscalar and axial- vector mesons appear as poles 
in the inhomogeneous Bethe-Salpeter amplitude associ- 
ated with the axial- vector vertex, . An exact form of 
the associated BSE is (q± = q ± P/2, etc.) 

rg(fc; P) = Z 2757ai -zJ g 2 D aP (k - q )^- la S f (q + ) 
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dq 
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x-A&(fc,g;P), 



(3) 



where Ag®^ is a four-point function, completely defined 
H via the quark self-energy and hence the quark-gluon 
vertex, r M . Crucially, Ag?g satisfies a Ward-Takahashi 
identity, whose solution provides a symmetry-preserving 
Ansatz consistent with T M . Bound-state masses in the 
++ , 1 __ , l" 1 channels arc, respectively, obtained by 
replacing Z 2 r y^ tl by the structures Id, 7/*, 75^, each 
multiplied by an appropriate renormalization constant. 



A prediction for the spectrum therefore follows once 
the gap equation's kernel is specified and the Ward- 
Identity solved for A{^. The kernel may be rendered 
tractable by writing fl7l[il| 

Z l9 2 D pa {t)T a {q, q + t)= Git 2 ) D i ™ e (t)Z 2 r a (q, q + t), 

(4) 

wherein D p ™ e is the Landau-gauge free-gauge-boson pro- 
pagator, Q is an interaction model and T a is an Ansatz for 
the quark-gluon vertex. For the interaction, we employ 



G(s) Sir 2 
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where: 7m = 12/25, A QCD = 0.234 GeV; r = e 2 - 1; 
and F{s) = {1 - cxp(-s/[4m 2 ])}/s, m t = 0.5 GeV. This 
form preserves the one-loop renormalization-group be- 
havior of QCD in the gap- and Bethc-Salpeter-equations, 
and is consonant with modern DSE- and latticc-QCD 
results. Detailed explanations of its development and 
capability as a tool in hadron physics are presented 
in Refs. [9j, Il9i. We employ the renormalization proce- 
dures of Ref. [17j and the same renormalization point, 
C = 19 GeV. N.B. This study is the first to employ the 
fully-renormalized forms of the vertex-dressed gap equa- 
tion and the symmetry-preserving Bethe-Salpeter equa- 
tions that are derived from it. 

We adapt the vertex explained in Refs. [Hl[2(|[2l]], vrz - : 



f M (p! , P2 ) = r^ c ( Pl , P2 ) + r* cm (pi , P2 ) 



(6) 



K° (Pi , P2 ) = *S A [p\ , p\ ) 7/1 + 2i p . [i 7 ■£ A A (p 2 , p\ ) 

+&b(pIpD] , (7) 

where ^{p 2 ,p 2 ) = [<i>(p 2 ) + 0(p§)]/2, A (p 2 ,p 2 ) = 
[4>{p\) — ~P 2 ]i = Pi +P2', and the anomalous 

chromomagnctic moment piece is 



r^ cm (Pi,p 2 ) = r 



acm i 



(Pl,P2) 



(P1.P2), (8) 



with (k = pi -p 2 , = Sft,, - k^kv/k 2 , a T p := T^a v ) 
rr n4 = [<?7 ■ * + hJ<r vp t»k p )r*(pi,P2) , (9) 

r^ cm5 = a^h-nton,!*) , (io) 



2-r 5 (pi,p 2 ) 

M(p 2 ,p 2 ) 



(11) 



T5 



nA B {p 2 ,p 2 ) and M(x,y) = [x + M(x) 2 
M(y) 2 }/(2[M(x) + M(y)}). 

Our Ansatz combines information and results from 
perturbative QCD, Dyson-Schwinger equation (DSE) 
and lattice-QCD studies. The structure and importance 
of the BC component has long been recognized (22|: 
its presence and momentum dependence are confirmed 
in comparisons between DSE- and lattice-QCD studies 
j23Tj25| . The ACM piece is novel, although its presence 
should long ago have been appreciated, given the impor- 
tance of the DCSB-induced A^-term in the BC compo- 
nent and the ACM term's identical origin. 
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Our introduction of the ACM term comes at a criti- 
cal juncture, since the developments in Ref. Q only now 
enable an exploration of its consequences for hadron ob- 
servables. Indeed, notwithstanding the large current- 
quark mass employed (115 MeV), a comparison between 
lattice-QCD estimates of T5(pi,j52) [HI and our ex- 
pression confirms the character of our Ansatz. For ex- 
ample, there is a two orders-of-magnitudc enhancement 
over the perturbative result, emphasizing the connection 
with DCSB, and semi-quantitative agreement between 
the momentum-dependence of the lattice result and our 
form. It is notable that in perturbation theory, whilst the 
on-shell anomalous chromomagnetic moment is always 
negative for a dressed-quark, the sign of the 75(^1,^2) 
contribution is positive in Landau gauge but negative in 
Feynman gauge (26|. To complete the picture, a further 
consideration of perturbation theory shows the vertex 
must also contain the T4-term because only then can the 
full vertex Ansatz reproduce the one- loop result fl3l.[26j. 
Moreover, Eq. (fTTj) is precisely the result one obtains on- 
shell in one-loop Landau-gauge perturbation theory. 

In the 0~ and 1 + channels the Ward-Takahashi iden- 
tity for the Bethe-Salpeter kernel is solved by 
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FIG. 1. Dashed curve ~ Dressed-quark anomalous chromo- 
magnetic moment distribution computed using 77 = 0.65; 
and solid curve - Electromagnetic moment distribution com- 
puted from the inhomogeneous Bethe-Salpeter equation for 
the dressed-quark-photon vertex. In the chiral limit and 
the absence of DCSB, both curves are zero. Moreover, ab- 
sent DCSB-induced terms in the dressed-quark-gluon vertex, 
Eq. ([6|, both curves would be an order-of- magnitude smaller 
. (The Euclidean constituent-quark mass M E := {p I p > 
0,p 2 = M 2 (p 2 )} = 0.35 GeV.) 



2A 5(S(Al) = \Tp(q + ,k+) +7 5 r /3 ( (? _,fc_)75] 

X 5-i(fc + )+Vi(-fc_) r ^ (fc;P) 

+ r*M(*P) s -i { _ q+) + s -i {q _ ) 

x fofcf/j (q+ , fc + ) 7 5 + f p (g_ , k- )] . (12) 

Given the vertex in Eq. ([5]) , we have now completely spec- 
ified equations in the 0~ and 1 + channels from which one 
may obtain bound-state masses and amplitudes. N.B. 
This solution of the Ward-Takahashi identity for the 
Bethe-Salpeter kernel is far more general than that pre- 
sented in Ref. Q. Kernels of equal simplicity and power 
for other channels are readily constructed by analogy. 



III. CHROMO- AND ELECTRO-MAGNETIC 
MOMENTS 



Before reporting results for the spectrum it is valuable 
to illustrate the novel DCSB content of Eq. ©. This may 
readily be accomplished via a single curve that depicts 
the dressed-quark anomalous magnetic moment distribu- 
tion, which is identically zero in the chiral limit fl3| . One 
characterizes dressed-quarks through a magnetic moment 
distribution because a confined quark does not possess a 
mass-shell 0, HH and hence one cannot unambiguously 
assign a single value to its anomalous magnetic moment. 

The magnetic moment distribution is computed as fol- 
lows. At each value of p 2 , we define spinors to satisfy a 
Dirac equation with the constant mass replaced by the 
running mass: m — > M(p 2 ) =: <;, and use the associated 



Gordon identity to write 

= S( P/ )[^ 1 (fc 2 ) 7Al + 1 a^KF^k^uipi). (13) 

Now, from Eqs. ^ - (|11[) , one finds an anomalous chro- 
momagnetic moment distribution 



'(0 



2^(3,7-1)+^) 



2^ 



(14) 



where a A 



S A (?,?)i 5 a = A a(<m0, etc. In Fig.rjjwc 
depict the result obtain using the interaction in Eq. . 

One computes the dressed-quark anomalous electro- 
magnetic moment distribution from the solution of the 
inhomogeneous Bethe-Salpeter equation in the color- 
singlet vector vertex. Details may be found in Ref. [l3| . 
which also explains the impact of DCSB in the dressed- 
quark-gluon vertex. The result is depicted in Fig.[TJ The 
marked similarity between this computed result and that 
produced by the Ansatz in Ref. plj is noteworthy. 



IV. MESON SPECTRUM 

We compute masses using the method detailed in 
Ref. [23] , which ensures one need only solve the gap and 
Bethe-Salpeter equations at spacelike momenta: a signif- 
icant numerical simplification. To explain, the inhomoge- 
neous BSE is solved for the complete amplitude in a par- 
ticular channel on a domain of spacelike total-momenta, 
P 2 > 0. Any bound-state in that channel appears as 



a pole in the solution at P 



Denoting the 
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FIG. 2. Illustration of the procedure used to determine meson 
masses. Solid curve - ai-meson, nonperturbative kernel; dot- 
dash-dash - a\, kernel derived from Eq. © only (Ball-Chiu, 
BC); and dash - ai, kernel derived from just the first term in 
Eq. (0 (1BC, a minimal renormalization improvement pl| of 
the leading-order - RL, rainbow-ladder - kernel [T3])- Dot- 
dash curve - p-meson, nonperturbative kernel; Dot-dash-dot 
- p, BC-kernel; and dotted - p, lBC-kernel. Points - values of 
1/T(k = 0;P 2 ) in the given channel computed with the kernel 
described. Pade approximants are constructed in each case; 
and the location of the zero is identified with (-mmesou)- 



leading Chebyshev moment of the amplitude's dominant 
Dirac structure by T(fc;P), then l/T(k = 0;P 2 ) exhibits 
a zero at (— ra^ cson ). The location of that zero is deter- 
mined via extrapolation of a Pade approximant to the 
spacelike-bchavior of l/r(fc = 0;P 2 ). This is illustrated 
for the p- and ai-channels in Fig. [5] 

Our results are listed in Table U wherein the level of 
agreement between Cols. 3 and 4 illustrates the efficacy 
of the method we're using to compute masses: no dif- 
ference is greater than 1%. Next consider m a and com- 
pare Cols. 1-3. It is an algebraic result that in the RL- 
truncation of QCD's Dyson-Schwingcr equations (DSEs), 
m a « 2AI , where M is a constituent-like quark mass [28j [. 
On the other hand, incorporating the quark mass func- 
tion into the Bethe-Salpeter kernel via r^ c generates a 
strong spin-orbit interaction, which significantly boosts 
rria [8[ . This feature is evidently unaffected by the inclu- 
sion of r^ cm ; i.e., those terms associated with a dressed- 
quark anomalous chromomagnctic moment. Since we de- 
liberately omit terms associated with pion final-state in- 
teractions in our nonperturbative kernel, it is noteworthy 
that m a in Col. 1 matches estimates for the mass of the 
dressed-quark-core component of the er-meson obtained 
using unitarized chiral perturbation theory (29l l30j . 

Now compare the entries in Rows 2, 4-6. The p- and 
ai-mesons have been known for more than thirty years 
and are typically judged to be parity-partners; i.e., they 
would be degenerate if chiral symmetry were manifest 
in QCD. Plainly, they are not, being split by roughly 
450 MeV (i.e., > m p /2). It is suspected that this large 
splitting owes to DCSB: hitherto, however, no symmetry- 
preserving bound-state treatment could explain it. This 
is illustrated by Cols. 3, 4, which show that whilst a 



TABLE I. Col. 1: Spectrum obtained with the full nonpertur- 
bative Bethe-Salpeter kernels described herein, which express 
effects of DCSB: u = 0.5 GeV; Duj = (0.52 GeV) 3 ; and we 
assume isospin symmetry, with rnj — = m = 3.7 MeV, 
£ = 19 GeV. The method of Ref. [27J was used: the error re- 
veals the sensitivity to varying the order of Pade approximant. 
We compute f n = 0.091 GeV using Eq. (15) in Ref. @. Col. 2 
- Experimental values; computed, except m CT , from isospin 
mass-squared averages [3lJ • Col. 3 - Masses determined from 
the inhomogeneous BSE at leading-order in the DSE trunca- 
tion scheme of Ref. [T3| using the interaction in Ref. || (with 
this simple kernel, the Pade error is negligible); and Col. 4 - 
results in Ref. [33], obtained directly from the homogeneous 
BSE at the same order of truncation. In the last two cases, 
m = 5 MeV, u = 0.5 GeV, Dlj = (0.794 GeV) 3 . 
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RL-Pade 


RL-direct 




0.138 


0.138 


0.138 


0.137 


m p 


0.84 ± 0.03 


0.777 


0.754 


0.758 


m a 


1.13 ±0.01 


0.4 - 1.2 


0.645 


0.645 


m ai 


1.28 ±0.01 


1.24 ±0.04 


0.938 


0.927 


m bl 


1.24 ±0.10 


1.21 ±0.02 


0.904 


0.912 


m ai - m p 


0.44 ± 0.04 


0.46 ± 0.04 


0.18 


0.17 


m bl — m p 


0.40 ±0.14 


0.43 ± 0.02 


0.15 


0.15 



good estimate of m p is readily obtained at leading-order 
in the systematic DSE truncation scheme of Ref. [l0| , the 
axial- vector masses are much underestimated. The flaw 
persists at next-to-leading-order [ll], [13] ■ 

Our analysis points to a remedy for this longstand- 
ing failure. Using the Poincare-covariant, symmetry pre- 
serving formulation of the meson bound-state problem 
enabled by Ref. @, with nonperturbative kernels for the 
gap and Bethe-Salpeter equations, which incorporate and 
express effects of DCSB that are impossible to capture 
in any step-by-step procedure for improving upon the 
rainbow-ladder truncation, we provide realistic estimates 
of axial- vector meson masses. In obtaining these results 
we found that the vertex Ansatz used most widely in 
studies of DCSB, , is inadequate as a tool in hadron 
physics. Used alone, it increases both m p and m ai but 
yields m ai — m p = 0.22 GeV, qualitatively unchanged 
from the rainbow-laddcr-like result (see Fig.[T]). A good 
description of axial- vector mesons is achieved by includ- 
ing interactions derived from r^ cm ; i.e., connected with 
the dressed-quark anomalous chromomagnetic moment 
[l3l ]. Moreover, used alone, neither term in r^ cm can 
produce a satisfactory result. The full vertex Ansatz and 
the associated gap and Bethe-Salpeter kernels described 
herein are the minimum required. 

In preparing Row 5 we obtained additional informa- 
tion. The leading-covariant in the 61-meson channel 
is 75^. The appearance of k p suggests that dressed- 
quark orbital angular momentum will play a signifi- 
cant role in this meson's structure, even more so than 
in the ai-channel for which the dominant covariant is 
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757m- (NB. In a simple quark-model, constituent spins 
are parallel within the a\ but antiparallel within the 
61. Constituents of the b\ can therefore become closer, 
so that spin-orbit repulsion can exert a greater influ- 
ence.) This expectation is realized in the result that 
mb 1 is far more sensitive to the interaction's range pa- 
rameter, w, than any other state, increasing rapidly with 
decreasing uj. Such behavior is readily understood. The 
kernel's nonpcrturbativcly-induced spin-orbit interaction 
acts over a length-scale characterized by r c := 1/ui. As r c 
is reduced with D fixed, the effective range of spin-orbit 
repulsion is also reduced, and the mass therefore drops. 

V. EPILOGUE 

Our results rest on an Ansatz for the quark-gluon ver- 
tex. The best available information was used in its con- 
struction. Improvement is nonetheless possible, which 
will involve elucidating the role of Dirac covariants not 
yet considered and of resonant contributions; viz., meson 
loop effects that give widths to some of the states we've 
considered. In cases for which empirical width-to-mass 
ratios are < 25%, we judge that such contributions can 
reliably be obtained via bound-state perturbation theory 
[33| . Contemporary studies indicate that these effects re- 
duce bound-state masses but the reduction can uniformly 



be compensated by a modest inflation of the interaction's 
mass-scale [28|,[34j|, so that the masses in Tableware semi- 
quantitativcly unchanged. The case of the cr-meson is 
more complicated. However, we predict a large mass for 
this meson's dressed-quark core, which leaves sufficient 
room for a strong reduction by resonant contributions 
[29I . l3(il ] . Balance requires we note that one cannot yet 
completely exclude the possibility of strong couplings be- 
tween other channels, consequences of which for under- 
standing the spectrum are illustrated in Ref. [3oT ] . 

We furnished a continuum framework for computing 
and explaining the meson spectrum. It combines a vera- 
cious description of pion properties with estimates for 
masses of light-quark mesons heavier than m p . Our 
method therefore offers the promise of a first reliable 
Poincare-invariant, symmetry-preserving computation of 
the spectrum of light-quark hybrids and exotics; i.e., 
those putative states which are impossible to construct 
in a quantum mechanics based upon constituent-quark 
degrees-of-freedom. So long as the promise is promptly 
fulfilled, the approach will provide predictions to guide 
the forthcoming generation of facilities. 
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